Atmospheric leptons 
the search for a prompt component 
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The flux of high-energy (>GeV) neutrinos consists primarily of those produced by cosmic-ray 
interactions in the atmosphere. The contribution from extraterrestrial sources is still unknown. 
Current limits suggest that the observed spectrum is dominated by atmospheric neutrinos up to at 
least 100 TeV. The contribution of charmed hadrons to the flux of atmospheric neutrinos is important 
in the context of the search for astrophysical neutrinos because the spectrum of such "prompt" 
neutrinos is harder than that of "conventional" neutrinos from decay of pions and kaons. The prompt 
component therefore becomes increasingly important as energy increases. This paper reviews the 
status of the search for prompt muons and neutrinos with emphasis on the complementary aspects 
of muons, electron neutrinos and muon neutrinos. 



I. INTRODUCTION 

Before experimental discovery of charmed hadrons at 
accelerators in the mid-1970s there was intense interest in 
using atmospheric muons to find evidence for production 
of heavy, short-lived hadrons. For example, a highlight 
of the 1973 cosmic-ray conference at Denver was an up- 
date of measurements made over several years with the 
underground muon spectrometer at Park City, Utah [![. 
At first, the observed angular dependence of multi-TeV 
muons had appeared to be more isotropic than could be 
explained solely by production through the pion and kaon 
channels, which are strongly enhanced at large zenith an- 
gles. With an improved understanding of the overburden 
and the detector response, however, it was finally con- 
cluded that the Park City data were consistent a "con- 
ventional" origin from decay of pions and kaons. An 
isotropic "prompt" component was not manifest for en- 
ergies below 10 TeV. 

Production of charmed hadrons has now been mea- 
sured over a large range of energy at accelerators. The 
production cross section increases significantly from ap- 
proximately ~ 1 /A> at y/s w 10 GeV to several mb at 
TeV 0. There is still not full coverage of phase 
space for charm production, however. In particular, the 
level of "intrinsic" charm Q production is still uncer- 
tain. The SELEX measurement shows a large asym- 
metry in the ratio of charm to anti-charm baryons pro- 
duced by baryon beams on a fixed target, while little or 
no asymmetry is observed in a pion beam. This obser- 
vation indicates some level of intrinsic charm in which 
the valence quarks of the projectile pick up a charmed 
quark. Charmed hadrons produced as fragments of the 
incident nucleon beam will contribute disproportionately 
to the spectrum of atmospheric leptons because of the 
steep cosmic-ray energy spectrum. Thus, even if intrin- 
sic charm contributes less to the total cross section for 
producing charm than production via QCD processes, it 
may have a significant effect on the prompt contribution 
to atmospheric muons and neutrinos. 

In addition to the intrinsic interest in identifying the 



charm contribution to the fluxes of atmospheric muons 
and neutrinos, there is another, perhaps more important, 
reason for trying to measure this component. That is be- 
cause of the relevance of prompt neutrinos in the search 
for neutrinos of astrophysical origin. Like a flux from un- 
resolved extra-galactic neutrino sources, the prompt con- 
tribution is isotropic for E v < 10 7 GeV. It is also harder 
by one power of energy than the spectrum of conventional 
atmospheric neutrinos. For these reasons, prompt neu- 
trinos constitute an important background for neutrino 
astronomy. 

The paper begins in § [IT] with a discussion of the ingre- 
dients needed to calculate fluxes of atmospheric muons 
and neutrinos, including relevant analytic approxima- 
tions and the primary cosmic-ray spectrum. Section Mil 
reviews models for charm production and corresponding 
predictions for fluxes of muons and neutrinos. In § IIVI we 
calculate the fluxes of conventional atmospheric muons 
and neutrinos and compare them with the charm contri- 
bution. The effect of the knee of the primary spectrum 
is included. The predictions are illustrated in § [V] with 
approximate calculations of the event rates for detector 
with a gigaton target volume like IceCube 0-0|- The 
concluding Section I VII comments on the current status 
and prospects for detection of prompt leptons in the near 
future. 



II. ATMOSPHERIC MUONS AND NEUTRINOS 

The two main ingredients in the calculation of at- 
mospheric neutrinos are the primary spectrum and the 
hadronic physics of meson production in hadronic inter- 
actions. Because production of pions, kaons and charmed 
hadrons occurs at the nucleon level, what is most relevant 
is the primary spectrum of nucleons per GeV/nucleon. 
Composition comes in through the ratio of protons to to- 
tal nucleons, which determines the charge ratio of muons 
and particle/anti-article ratio for neutrinos. 
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FIG. 1: Left: All-particle spectrum from Ref. || where references to the data are given. Right: Spectrum of nucleons for 
several assumptions (see text for explanation). 



A. Primary spectrum 



spectrum is given by 



For illustration I use a phenomenological model of the 
primary spectrum with three populations of particles and 
five nuclear components @, as shown in Fig. [1] (Left). 
There are two basic assumptions. First, it is assumed 
that all energy dependence (whether from acceleration or 
propagation) depends only on how particles are affected 
by their magnetic environment. As a consequence, each 
nuclear component (mass number Z and total momen- 
tum per particle P) depends on magnetic rigidity (R) 
in the same way, where R = Pc/Ze. Peters Q pointed 
out the consequence of this assumption for the primary 
composition in the region of the knee of the spectrum, 
namely, that, when expressed in terms of total energy 
per particle, protons would steepen first followed by he- 
lium and then by nuclei with successively higher charge. 

The other assumption, following Hillas [lfj, is that 
three populations of particles are sufficient to charac- 
terize the entire cosmic-ray spectrum. This is almost 
certainly an oversimplification. A more realistic picture 
would likely involve many individual sources injecting 
particles at various distances and times, as in the model 
for galactic cosmic rays of Blasi and Amato [11]. Thus 
the three populations represent three classes of sources: 

1. Particles accelerated by supernova remnants in the 
galaxy, 

2. A higher energy galactic component of uncertain 
origin, and 

3. Particles accelerated at extra-galactic sources. 
The contribution of nuclei of mass A4 to the all-particle 



ME) 



EdNi 
dE 



E- 



x exp 



E 



ZiR c 



(1) 

where E is the total energy per nucleus. 

The spectral indices for each group and the normal- 
izations are given explicitly in Table |T] The parameters 
for Population 1 are based on fits to spectra of nuclear 
groups measured by CREAM [l2|, [H| , which we assume 
can be extrapolated to a rigidity of 4 PV to describe the 
knee. This is an unverified simplifying assumption that 
needs to be checked by measurements in the PeV range. 
In Eq. Q] fa is dN/ din E and 7, is the integral spectral 
index. The subscript i = 1, 5 runs over the standard five 
groups (p, He, CNO, Mg-Si and Fe), and the all-particle 
spectrum is the sum of the five. 

The spectrum of nucleons as a function of energy per 
nucleon corresponding to Eq. [T] is given by 
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Because of the steep cosmic-ray spectrum, protons are 
relatively more important and heavy nuclei less impor- 
tant in the spectrum of nucleons (as a function of En — 
E ioi /A) than in the all particle spectrum. 

The spectrum of nucleons is plotted for several assump- 
tions in Fig. [1] (right). The straight solid line shows a 
simple E~ 2 - 7 spectrum of nucleons to guide the eye. The 
straight dotted line shows the spectrum below 100 TeV 
recommended in 2001 [l4{ as a standard for use in calcu- 
lating fluxes of atmospheric leptons up to 10 TeV. At low 
energy the fit was based on measurements of BESS (lijij 
and AMS [16]. The spectral index used for protons at 
the time was rather steep (2.74), based on the measure- 
ments of BESS and AMS below 200 GeV. Recent results 
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TABLE I: Cutoffs, integral spectral indices and normalizations constants atj for Eq. [T] 



of PAMELA [17j show that the spectrum of protons hard- 
ens above 200 GeV. Two options were given for helium, 
which contributes of order 25% of the spectrum of nucle- 
ons. The high option for helium (with an integral spec- 
tral index of 7 = 1.64) suggested by emulsion chamber 
measurements in the 10 TeV range at the time fill Il9l ] 
has since been confirmed by ATIC 0, CREAM [l2| and 
PAMELA I7i]. A version of the spectrum of Ref. GHLS 
is used in the standard calculations of the flux of atmo- 
spheric neutrinos by Honda et al. [2l| and by the Bartol 
group [22| . The spectrum of Honda et al. (as described 
in [23j) uses a harder spectrum for hydrogen (2.71 in- 
stead of 2.74) above 100 GeV. Their overall spectrum is 
nearly constant at 7 = 1.69 from 200 GeV to 50 TeV with 
a fraction of helium that increases from 20% to 25% in 
the same region. The spectral index of the spectrum of 
nucleons in the model of Ref. Q is nearly constant at 
7 = 1.63 from 200 GeV to 50 TeV with a corresponding 
increase in the contribution of helium from 22% to 30%. 
The contribution from heavier nuclei is at the level of 
10%. 



Particle (a): \ n ± K± K° L Charm 



The other lines in the spectrum of nucleons all include 
the effect of the knee in the cosmic-ray spectrum in one 
way or another. The heavy solid curve is the nucleon 
spectrum corresponding to the model in Table Q] The 
nearby pink dash-dot is an analytic approximation to 
that model, which is described below in Eq. [T3] The 
strong knee around 1 PeV is the consequence of the in- 
creasing fraction of heavy nuclei in the model. In addition 
to the model of Ref [81, Fig. Q] also shows the polygonato 
model [24[ without any contribution from nuclei heavier 
than iron. Each nuclear component in the model steep- 
ens by 8 = 1.9 at a rigidity of 4.49 PV. The effect of the 
knee begins to show up in the nucleon spectrum already 
somewhat below one PeV. Using the rule of thumb that 
there is on average a factor of ten between the parent 
cosmic ray energy and the secondary leptons, taking ac- 
count of the steepening of the spectrum will be important 
for muon and neutrino energies of 100 TeV and above, 
which we is discussed in § [TV] The double dotted line 
that steepens from a differential index of —2.7 to —3.0 at 
5 x 10 6 GeV is the primary spectrum used in the charm 
calculation of Ref. TIG. 



B. 



e a (GeV): 115 850 205 ~ 3 x 10' 
TABLE II: Characteristic energies. 



Hadron production and decay kinematics in the 
atmosphere 



The phenomenology of atmospheric leptons depends 
on the production of pions, kaons and heavier hadrons 
by interactions of cosmic-rays in the atmosphere and on 
the kinematics for the relevant decay channels. Produc- 
tion and subsequent decay occur through generation by 
a steep spectrum of primary and secondary cosmic rays 
in the atmosphere. The competition between reinterac- 
tion and decay of unstable hadrons depends on density 
and altitude. In the framework of a set of analytic ap- 
proximations for solution of the cascade equations, the 
essential dependence on energy and zenith angle comes 
through the critical energy defined as 
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where the index 1 indicates the hadron (7r ± , K^, Kl 
or charmed hadron), r, is the meson lifetime, h n is the 
scale height of the atmosphere and 8* is the zenith angle 
(*corrected for curvature of the Earth for 9 > 70°). Val- 
ues of the important characteristic energies are given in 
Table EH 

For a power-law spectrum of primary nucleons, the ex- 
pression for the lepton spectrum factorizes into a prod- 
uct of the primary spectrum and an expression that re- 
flects the properties of production of secondary hadrons 
by the cosmic-ray spectrum and their subsequent decay 
to muons and neutrinos. 
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The ^4-factors in Eq. [4] are a product of the spectrum 
weighted moments for production of mesons by nucleons 
times the spectrum weighed moments of the meson decay 
distributions, which include both the decay kinematics 
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and the branching ratios. For a power-law spectrum of 
decaying pions with a differential spectral index a the 
decay factor is 
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for muons and neutrinos respectively. In the low energy 
limit, the spectral index is the same as that of the pri- 
mary spectrum of nucleons, but at high energy the spec- 
trum of the decaying pions is one power steeper because 
of the competition between decay and reinteraction. Low 
and high are with respect to the critical energy e„/ cos(#). 
The ratio r n — m„/m^ = 0.5731. 

The forms for two-body decay of charged kaons are 
the same, but the mass ratio factor is much smaller: 
tk = 0.0458. The larger critical energy for charged kaons 
leads to an increase in the contribution of kaons with 
increasing energy for both muons and neutrinos. The 
differences between the kinematic factors for two-body 
decay to neutrinos and muons amplifies the importance 
of the kaon channel for neutrinos. At high energy, in the 
TeV range and above, charged kaons account for about 
80% of muon neutrinos as compared to 25% of muons. 

Each term in Eq. |4] is a form that combines the low 
energy and high energy solutions to the cascade equa- 
tion respectively for pions, kaons and charmed hadrons 
in the atmosphere. The numerator is a product of the 
spectrum weighted moment for meson production and 
the spectrum weighted moment of the decay distribution 
to [25| with a — 7 + 1, where 7 is the integral spectral 
index of the spectrum of primary nucleons. The denom- 
inator governs the transition between the low and the 
high energy regimes. The forms for muons are similar. 
For low energy, meson decay dominates over reinteraction 
and the resulting lepton spectrum has the same shape 
as the primary spectrum of nucleons. Charmed hadrons 
are in the low-energy regime for -Ei op ton < 10 7 GeV. For 
high energy (-Ei ep ton > e Q /cos0), reinteraction of the 
hadron is more likely and the lepton spectrum becomes 
one power of energy steeper than the primary spectrum. 

In the high energy limit, the spectrum weighted mo- 
ment for meson decay has to be evaluated on the steeper 
spectrum, and the attenuation lengths for reinteraction 
come into play. The By quantities in the denominators 
are the product of the ratios of low-energy to high en- 
ergy decay distributions combined with the function of 
attenuation lengths that accounts for cascading of the 
mesons [26j |. Explicitly, for neutrinos 
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The forms for kaons are the same as functions of tk and 
Ak- The dependence of 7 on energy in the case of a 



non-power law primary spectrum needs to be accounted 
for. 

For a power-law primary spectrum of nucleons and as- 
suming Feynman scaling for hadron production, the cas- 
cade equations can be solved analytically as in Eq. 2) 
The primary spectrum can always be described locally as 
a power law, and similarly the hadronic interactions can 
be written in terms of the scaled energy (x = E a / En) 
for a given primary energy per nucleon. In both cases the 
dependence on energy is sufficiently gradual that the ap- 
proximate analytical forms can be used for quantitative 
calculations if the slow variation with energy is accounted 
for. This approach is taken in the calculation of Thun- 
man, Ingelmann and Gondolo (TIG) [271, which I follow 
here. They define energy-dependent Z-factors as in the 
following example for nucleons producing charged pions: 



Z N7r ±(E) — 



(E') X N {E) dn„±(E',E) 



4> N {E) \ N {E>) 



dE 
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Here Ajy(B) is the nucleon interaction length and dn^ is 
the number of charged pions produced in dE by nucle- 
ons of energy E', and 4>n{E) is the spectrum of nucleons. 
The energy-dependent Z-factors are then used in applica- 
ble version of Eq. |4]to evaluate the lepton spectrum. This 
approximation is valid to the extent that the energy de- 
pendences are gradual. They showed that the numerical 
approximation based on the spectrum weighted moments 
taken from the interaction model used in their Monte 
Carlo produced similar results to the full Monte Carlo. 
The advantage is that the analytic approximations can 
be tuned to match well to a particular model of hadronic 
interactions. They can then be used to extend a Monte 
Carlo calculation to arbitrarily high energies without the 
statistical problems that arise from the fact that mesons 
usually interact rather than decay at hi gh e nergy. This 
method has been used recently in Ref. [28] to comple- 
ment Monte Carlo calculations of atmospheric leptons in 
the region of the knee of the cosmic-ray spectrum. 



III. EXPECTATIONS FOR PROMPT LEPTONS 

Limits on the prompt contribution to atmospheric 
muons are conventionally expressed in terms of the ratio 
R c of charm production to pion production by rewriting 
Eq. H as 



cos6En/e n 



A 



1 + B K/J , cosOE^/ex 



This form applies for lepton energy < 10 7 GeV where the 
charm contribution is isotropic. The current upper limit 
from LVD [29] is R c < 2xl0 -3 assuming a differential 
primary spectrum oc E~ 2 ' 77 . The primary spectrum of 
nucleons used in this analysis is less steep by about 0.1, so 
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FIG. 2: Predictions of three models for the flux of prompt 
muons. See text for discussion of the rescaled plots. 



the LVD limit at 10 TeV (for example) would be reduced 
by a factor of two for comparison with the models as 
discussed here. In what follows I will compare the current 
models for prompt leptons to R c ~ 10~ 3 considered as 
an experimental upper limit. 



A. Calculations of charm production 

Reference [27| uses a simple Monte Carlo to generate 
the distribution of hadronic interactions and decays in 
the atmosphere and Pythia [3(3] to generate the secondary 
hadrons at each interaction point. Charm production is 
calculated within Pythia using first order QCD matrix 
elements to calculate cc production by gluons and by 
quarks. A renormalization factor K — 2 is used to rep- 
resent higher order QCD effects. In addition to their 
Monte Carlo calculation, They parameterized their re- 
sults for the charm contribution in a form similar to the 
third term of Eq. 2] as 



cf>c(E) = N 
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= 1.77 below and 
for this model is 



2.0 

2 x 



with A w 3 x 10 7 GeV and 7 
above 10 6 GeV. The value of R 
10~ 5 at 10 TeV. 

The reference calculation used for evaluation of the 
atmospheric neutrino background in IceCube at high en- 
ergy is that of Enberg, Reno & Sarcevic (ERS) [31|. It 
is a QCD calculation that gives results somewhat higher 
than some previous calculations [27l [32| and lower than 
an earlier NLO-QCD calculation [33]. The ERS calcula- 
tion assumes the same primary spectrum of nucleons as 



in TIG [23], which is shown in the right panel of Fig. [T] 
The shape of the ERS calculation is similar to that of 
TIG, and its central value is approximately a factor of 
two higher. ERS assign an uncertainty range of approx- 
imately ±50% to their calculation. The value of R c for 
this model is ~ 10~ 4 at 10 TeV. 

The "Recombination quark-parton model" 
(RQPM) [34[ embodies the idea of intrinsic charm. 
The underlying concept is that-also for heavy quarks- 
there is a process of associated production in which the 
cc pair produced when the projectile proton fragments 
can rccombine with valence quarks (di-quarks) and 
with sea quarks to produce charmed hadrons, includ- 
ing charmed hyperons. For example, in this model 
the process p —> A+ + D° would be expected at a 
level (m s /m c ) 2 relative to associated production of 
strangeness, p — > AK + . The parameters of the model 
are adjusted to fit the then available data on charm 
production in Ref. [35], and a parameterization of the 
prompt muon flux in the energy range from 5 TeV to 
5000 TeV is given. The value of R c for this model 
is ~ 8 x 10~ 4 at 10 TeV. The RQPM model is close 
to the LVD upper limit, but still allowed by it. It is 
interesting that the recent IceCube limit on the prompt 
contribution to ^-induced upward muons [3(| is just at 
the level of the RQPM model. 

Figure O compares the predictions of these three mod- 
els for the prompt contribution to the atmospheric muon 
flux. An important point to note in comparing the curves 
in this figure is that the spectrum also depends on what 
is assumed for the primary spectrum. The TIG and ERS 
use the same assumption as each other, in which the knee 
is probably too high in energy. The primary spectrum 
used in Ref. [35| is taken from a model [13] in which the 
knee is attributed to energy losses in the sources (photo- 
disintegration for nuclei and photo-pion production for 
nuclei. In this case, the knee is a function of E/A rather 
than rigidity for the nuclei, and the knee occurs at higher 
energy per nucleon for protons than for helium because 
of the relatively high threshold for photo-pion production 
as compared to photo-disintegration. This spectrum also 
has a knee at rather high energy and, in addition, ap- 
pears to be anomalously high even in the few TeV range, 
as shown in Fig. [TJ The broken lines labeled "rescaled" 
in the figure are estimates of what the ERS and RQMP 
models would give for the prompt muon flux if the spec- 
trum of Ref. [8[ had been used. This estimate is obtained 
by multiplying by the ratio at = 10 x E v of the nu- 
cleon flux used here [|| to those used for ERS [3l| and 
RQPM [H]. 



IV. ATMOSPHERIC LEPTONS INCLUDING 
THE KNEE 



My presentation at ISVHECRI-2012 showed a series 
of figures in which the contribution of charm was shown 
separately from the "conventional" atmospheric leptons 
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FIG. 3: Muon spectra including prompt contribution. Left: prompt component from ERS model (rescaled); Right: prompt 
component from RQPM model (rescaled). 



from decay of pions and kaons. In those figures, a sim- 
ple power law spectrum was assumed for evaluating the 
fluxes of the conventional leptons, and the charm contri- 
butions were taken directly from the models. Here I take 
account of the knee in the spectrum, and I make a prelim- 
inary effort to present a consistent representation of the 
prompt component by rescaling their fluxes to the same 
spectrum of nucleons used for the conventional leptons, 
as described in the previous paragraph. 

Since the primary spectrum is no longer a power law, 
the spectrum-weighted moments depend on energy. To 
evaluate energy-dependent Z-factors, Eq. [8] is used. For 
simplicity here only the energy dependence of the spec- 
trum is considered, and standard values of nucleon in- 
teraction and attenuation lengths are used. In addition, 
a scaling form for meson production is assumed. The 
goal here is to demonstrate the effect of the knee in the 
cosmic-ray spectrum on the lepton fluxes. 



A. Approximation for pion and kaon production 



B. Approximation for nucleon spectrum 

The nucleon spectrum with the knee as parameterized 
in Table U and Eq. [2] can be approximated well (better 
than 10% to 30 PeV) with the standard two-power-law 
form of Ref. [3i| to describe the knee, a form which is 
also used in Ref. [24j]. Specifically, 

E^- = const x £TT (1 + {E/E*Y)- 5 '\ (13) 
aE 

with 7 = 1.64, 5 = 0.67, E* = 9.E + 5 GeV and 
e = 3.0. This approximation locates the knee in the nu- 
cleon spectrum just below a PeV, and the slope steepens 
from an integral spectral index of 1.64 below the knee 
to 2.31 after the knee. The normalization constant is 
10290. m -2 sr _1 s -1 . The steepening in the nucleon spec- 
trum is a consequence of the steepening of the all-particle 
spectrum amplified by the increasing fraction of nuclei in 
the all-particle spectrum. Equations [T^l H2] and [T3] are 
combined and integrated numerically to obtain energy- 
dependent Z-factors. 



Explicit approximations for scale-independent meson 
production forms are given in Ref. [38[ . The approxima- 
tion for nucleons to produce charged pions is 



Fnk{E v l 'En) — E n 



An 1T (E 1T /E N ) 
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with a similar form for production of kaons. With these 
scaling forms for particle production, the integral in Eq.[5] 
can be rewritten as 



r 1 dx Mg/g) 
Znv± - l 4> N (E) Fn ^ 



(12) 



C. Atmospheric leptons including the knee 



The fluxes of /i ± , 



from decay of pions and 



kaons are obtained using the energy-dependent Z-factors 
to evaluate Eq. [4] The energy-dependence of the spectral 
index in the meson decay factors that appear in Eq.Ulare 
also accounted for by using the local (energy-dependent) 
integral spectral index of the nucleon spectrum, which 
steepens gradually from 7 = 1.64 to 2.31 through the 
knee region. 

Calculation of the flux of electron neutrinos requires 
tracing the contributions of Ke3 decays of both charged 
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FIG. 4: Neutrino spectra including the prompt contribution. Left: -\- ; Right: v e + v a . 



and neutral kaons. This has been done by taking account 
of the neutron/proton ratio of the primary nucleons and 
tracking separately the production of KK pairs and pro- 
duction of kaons in association with A and £ hyperons. 
An approximate value of n/p — 0.54 has been used Q. 
As noted in the previous section, associated production 
of kaons by dissociation of an incident nucleon into a kaon 
and a hyperon is a prototype for intrinsic charm. 

Results for the lepton fluxes are shown in Figs. [3] and 2] 
The vertical muon fluxes are shown in Fig. [3l comparing 
the ERS model for prompt muons (left) to the RQPM 
model (right). Below the knee the RQPM prompt flux 
is rescaled down by about 15%, while the ERS prompt 
muon flux is rescaled up by about 20%. The rescaled 
ERS prompt muons are repeated on the right panel for 
comparison. In the region between 100 TeV and 1 PeV 
the rescaled RQPM prompt flux is a factor of 5 to 6 
higher than the rescaled ERS flux. The RQPM prompt 
component crosses the conventional flux around 300 TeV, 
as compared to a crossover at 2 PeV for ERS (left panel) . 

Figure [4] compares the situation for (left panel) with 
that for v e (right panel). The ERS model is shown for the 
prompt component. The RQPM flux is also repeated on 
both plots for reference. The conventional atmospheric 
v e flux is approximately a factor 20 lower than the con- 
ventional flux of f M , so the electron neutrino component 
is dominated by the prompt component at quite low en- 
ergy. 

V. EXPECTED RATES 

The fluxes described above can be used to estimate 
the rate of events in a kilometer scale detector. For at- 
mospheric muons the rate per year is simply the flux 
multiplied by 3 x 10 7 seconds/yr and by 10 10 cm 2 /km 2 . 



The corresponding integral rate of events I A ,(> E^) is 
shown in Fig. [5] 

The rate of neutrino-induced muons can be obtained 
in a similar way, with one additional step. It is necessary 
to calculate the effective area to convert the rate of neu- 
trinos with trajectories passing through the detector to a 
rate of neutrino-induced muons. Effective area is defined 
in such a way that 4>(E„,9) x A c s(E v ,8) is the rate of 
neutrino- induced muons per second per sr at zenith angle 
9. Explicitly 

A cS {E y ,6) = e(E thl 9)A(6)P„{E u ,E th ) (14) 
x exp{-*u(E v )N A X(0)}, 

where P(E V ) is the probability that a neutrino converts 
and produces a muon that reaches the detector with 
enough energy to be reconstructed. Absorption of neu- 
trinos in the Earth becomes significant in the 10 to 1000 
TeV range, first for vertically upward trajectories and 
for neutrinos with zenith angles ~ 120° around a PeV. 
An accurate calculation of A c g requires a detector sim- 
ulation. Here I use an estimate for an ideal km 2 de- 
tector from Ref. [4(| and estimate the rate of neutrino- 
induced muons in the zenith angle range from horizontal 
to —120°. The result is shown for the ERS assumption 
in the left panel of Fig. H>] 

Electron neutrinos must interact in the detector to 
be identified as cascades in the detector. Such cas- 
cades are virtually indistinguishable from neutral cur- 
rent interactions of muon or electron neutrinos of energy 
E v ~ E Vc jy 1 where y is the inelasticity of the neutral 
current neutrino interaction. The neutral current interac- 
tions of atmospheric make a comparable contribution 
to cascades for the conventional atmospheric neutrinos 
because the flux of is significantly higher than that of 
v e . For the charm component, however, the neutral cur- 
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FIG. 6: Estimate of the rate of atmospheric neutrino interactions per year in a km3 detector. Left: neutrino-induced muons 
per km 2 and with zenith angles from 90° to 120°; Right: electron neutrinos with vertices inside 1 km 3 . ERS is rescaled in both 
plots, as discussed in the text. 



rent contribution is relatively unimportant because of the 
equality of the fluxes of prompt z/ e and . The integral 
rate of v e interactions from all directions is 

R(> E v ) = 4nN x T x / a cc ^{E v )dE v , (15) 

where 4>{E U ) is the spectrum of v e + v e averaged over all 
directions and a cc is the charged current cross section, 
taken here from Ref. [4l|. N is the number of nucleons 
per cubic kilometer of ice and T=l year. The result- 
ing estimate of the rate of atmospheric v e interactions is 
shown in the right panel of Fig. [5] The plot includes the 



effect of neutrino shadowing by the Earth in the upward 
hemisphere [42j |. This amounts to a suppression for the 
whole sky of 9% at 100 TcV and 23% at a PeV. The plot 
shows about 5 electron neutrino interactions above 100 
TeV per km 3 of ice per year assuming ideal (full) effi- 
ciency. The spectrum is steeply falling so that less than 
one such event in ten years is expected above a PeV. 
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VI. SUMMARY COMMENTS 

In view of the recent observation by IceCube of two 
cascade-like events with observed energy just above 
a PeV [43[, together with the observation of several 
high=energy cascades reported at this meeting [44]], it is 
important to achieve a good understanding of the back- 
ground from atmospheric neutrinos in the energy region 
around 100 TeV and above. The work outlined in this 
paper is just a start. 

There are large uncertainties in the primary spectrum 
in the TeV range and above that need to be assessed. The 
model used here is just one possibility. A model for charm 
production that includes recent data from LHC [SH [46| 
in its fits needs to be developed. As an example, the 
framework for including charm within SIBYLL [48j ex- 
ists , but the parameters need to be tuned to data over 
a wide range of energy and phase space. An updated 



model could then be used for a full Monte Carlo sim- 
ulation including charm. The analytic approximations 
will remain an important tool to supplement the Monte 
Carlo, for example to track the consequences of uncer- 
tainties in the input spectrum and hadronic interactions 
for the expected fluxes [47[, as well as to parameterize 
and extrapolate the limited statistics of the full Monte 
Carlo [H. 
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